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1.  Foreword 

This  project  is  about  the  algorithm  development,  analysis,  implementation  and 
application  aspects  of  high  order  finite  difference  weighted  essentially  non-oscillatory 
(WENO)  schemes,  finite  volume  WENO  schemes,  discontinuous  Galerkin  finite  el¬ 
ement  methods  and  spectral  methods  for  solving  convection  dominated  problems 
requiring  long  time  integration  and  small  dissipation/dispersion  with  discontinuous 
or  high  gradient  solutions.  Algorithm  development  and  analysis,  investigation  about 
efficient  implementation  including  parallel  implementations,  and  applications  in  com¬ 
putational  fluid  dynamics,  computational  semiconductor  device  simulation  and  other 
areas,  are  performed.  The  achievement  strengthens  our  objective  to  obtain  power¬ 
ful  and  reliable  high  order  numerical  algorithms  and  use  them  to  solve  convection 
dominated  problems,  especially  those  of  army  interest. 


4.  Statement  of  the  Problem  Studied 

The  problems  studied  in  this  project  involve  numerical  solutions  of  convection 
dominated  partial  differential  equations.  These  problems  typically  have  solutions 
which  are  either  discontinuous,  or  with  discontinuous  derivatives,  or  containing  sharp 
gradient  regions  which  are  difficult  to  be  completely  resolved  on  today’s  computers. 
Our  objective  is  to  develop,  analyze  and  apply  numerical  methods  which  “capture” 
the  discontinuities  or  sharp  gradient  regions,  without  fully  resolving  them,  while  main¬ 
taining  nonlinearly  stable  transitions  for  these  discontinuities  or  sharp  gradients  and 
high  order  accuracy  in  the  smooth  part  of  the  domain.  High  order  accurate  finite 
difference  and  finite  volume  WENO  schemes,  finite  element  discontinuous  Galerkin 
methods,  and  spectral  methods  have  all  been  considered. 

Our  approach  is  to  explore  the  robustness  and  efficiency  of  high  order  numerical  al¬ 
gorithms  for  nonsmooth  problems  both  through  theoretical  guidance,  often  obtained 
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with  rigorous  proofs  on  simplified  model  problems,  and  through  numerical  experi¬ 
ments  on  real  application  problems.  We  do  not  try  to  modify  algorithms  just  for  the 
purpose  of  convergence  proofs,  if  such  modifications  are  not  justified  by  numerical 
experiments.  For  finite  difference  schemes,  we  are  exploring  the  very  efficient  WENO 
schemes  based  on  point  values,  numerical  fluxes,  and  nonlinearly  stable  high  order 
Runge-Kutta  time  discretizations.  For  finite  element  methods,  we  are  exploring  the 
Runge-Kutta  discontinuous  Galerkin  methods  of  Cockburn  and  Shu,  which  combine 
the  advantage  of  finite  elements  (weak  formulation,  automatic  energy  stability,  easy 
handling  of  complicated  geometry  and  boundary  conditions)  with  features  of  high 
resolution  finite  difference  schemes  (approximate  Riemann  solvers,  limiters).  Effec¬ 
tive  ways  to  handle  viscous  terms  and  higher  derivative  terms  are  being  investigated. 
For  spectral  methods,  we  are  exploring  reconstruction  techniques  of  Gottlieb  and  Shu 
to  apply  spectral  approximations  to  discontinuous  functions  and  still  obtain  uniform 
spectral  accuracy. 

We  have  been  continuing  on  the  study  of  efficient  and  high  order  finite  difference 
WENO  schemes  on  multiple  domains  with  overlaps,  which  will  be  useful  for  general 
problems  of  overlaying  domains  and  is  of  interest  to  Dr.  Rupak  Biswas  of  RIAGS 
and  Dr.  Roger  Strawn  of  US  Army  AFDD,  at  NASA  Ames  Research  Genter,  on 
the  investigation  of  developing  high  order  high  resolution  numerical  methods  for  the 
simulation  of  helicopter  rotor  blade  motion. 


5.  Summary  of  the  Most  Important  Results 

Research  has  been  performed  in  all  areas  listed  in  the  original  proposal,  and 
progress  and  results  consistent  with  the  original  objectives  have  been  obtained.  There 
are  54  publications  (among  them  32  appeared  in  refereed  journals,  6  appeared  in  con¬ 
ference  proceedings  and  book  chapters,  9  accepted  and  to  appear  in  refereed  journals, 
and  7  preprints  submitted  for  publications)  resulting  from  this  project,  see  Section  6 
for  a  list  of  them. 

S.  Gottlieb,  G.-W.  Shu  and  E.  Tadmor  have  reviewed  and  further  developed  a 
class  of  strong  stability  preserving  (SSP)  high  order  time  discretizations  for  semi¬ 
discrete  method  of  lines  approximations  of  partial  differential  equations,  [al]  (all 
the  numbering  of  references  are  according  to  that  of  Section  6).  Termed  TVD  (total 
variation  diminishing)  time  discretizations  before,  these  high  order  time  discretization 
methods  preserve  the  strong  stability  properties  of  first  order  Euler  time  stepping  and 
has  proved  very  useful  especially  in  solving  hyperbolic  partial  differential  equations. 
The  new  developments  include  the  construction  of  optimal  explicit  SSP  linear  Runge- 
Kutta  methods,  their  application  to  the  strong  stability  of  coercive  approximations, 
a  systematic  study  of  explicit  SSP  multi-step  methods  for  nonlinear  problems,  and 
the  study  of  the  strong  stability  preserving  property  of  implicit  Runge-Kutta  and 
multi-step  methods.  A  survey  of  the  SSP  time  discretizations  is  also  given  in  [bl]. 
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Discontinuous  Galerkin  (DG)  method  has  been  extensively  developed,  analyzed 
and  applied  during  this  period.  In  [a4],  Gockburn  and  Shu  have  given  an  exten¬ 
sive  review  of  the  state  of  the  art  of  this  method.  In  [a6,a7,a30,a31]  and  [dl2,dl4], 
Yan  and  Shu,  Levy,  Shu  and  Yan,  and  Xu  and  Shu  developed  and  analyzed  non- 
linearly  stable  local  discontinuous  Galerkin  (LDG)  methods  for  partial  differential 
equations  containing  third  and  higher  spatial  derivatives,  including  the  KdV  equa¬ 
tions,  time  dependent  bi-harmonic  equations,  K{m,n)  equations  with  compacton 
solutions,  KdV-Burgers  type  equations,  the  general  fifth-order  KdV  type  equations 
and  the  fully  nonlinear  iL(n,  n,  n)  equations,  the  Kuramoto-Sivashinsky  equations 
and  the  Ito-type  coupled  KdV  equations,  one  and  two  dimensional  generalized  non¬ 
linear  Schrddinger  equation  and  the  coupled  nonlinear  Schrbdinger  equation,  and  the 
two  dimensional  Kadomtsev-Petviashvili  equation  and  Zakharov-Kuznetsov  equation. 
These  LDG  methods  satisfy  cell  entropy  inequalities  and  are  nonlinearly  stable  in 
or  other  norms  for  quite  general  nonlinear  cases,  and  are  flexible  in  h-p  adaptivity  and 
efficient  for  parallel  implementation.  In  [al2,a25]  and  [d5],  Gockburn,  Luskin,  Shu 
and  Suli,  Ryan  and  Shu,  and  Ryan,  Shu  and  Atkins  have  explored  the  higher  order 
convergence  rates  in  negative  norms  for  discontinuous  Galerkin  methods  applied  to 
linear  hyperbolic  problems,  and  an  efficient,  local  post-processing  technique  which 
can  recover  (2A;-|-l)-th  order  of  accuracy,  instead  of  the  usual  (A;-|-l)-th  order,  in 
norm,  for  DG  solutions  to  linear  hyperbolic  problems  including  those  with  variable 
and  discontinuous  coefficients.  This  technique  can  also  recover  derivatives  of  the  nu¬ 
merical  solution  with  enhanced  order  of  accuracy.  Such  methods  are  expected  to  be 
extremely  useful  for  adaptive  computations  and  this  will  be  explored  in  the  future. 
In  [al8,d2],  Zhang  and  Shu  explored  a  method  to  analyze  the  convergence  and  rate 
of  convergence  for  discontinuous  Galerkin  methods  and  related  spectral  finite  volume 
method  (which  is  a  Petrov-Galerkin  method),  and  used  it  to  analyze  three  different 
formulations  of  the  discontinuous  Galerkin  method  for  solving  diffusion  problems,  as 
well  as  to  compare  discontinuous  Galerkin  method  with  the  spectral  finite  volume 
method.  In  [a23]  and  [d6,d7],  Qiu  and  Shu  developed  a  new  approach  using  weighted 
essentially  non-oscillatory  (WENO)  reconstructions  as  limiters  for  the  discontinuous 
Galerkin  methods  solving  hyperbolic  problems  containing  strong  discontinuities,  thus 
allowing  the  method  to  be  both  high  order  accurate  and  non-oscillatory  for  strong 
discontinuities.  A  Hermite  WENO  reconstruction  procedure,  which  relies  on  a  more 
compact  stencil  in  the  reconstruction  to  achieve  high  order  accuracy,  is  also  developed. 
In  [a28],  [b6]  and  [dl],  Gockburn,  Li  and  Shu,  and  Li  and  Shu  developed  the  locally 
divergence-free  discontinuous  Galerkin  method  for  solving  the  Maxwell  equations  and 
the  MHD  equations.  In  [a29],  Zhang  and  Shu  have  given  an  error  estimate  for  the 
fully  discrete  Runge-Kutta  discontinuous  Galerkin  method  applied  to  nonlinear  scalar 
hyperbolic  conservation  laws  with  smooth  solutions.  In  [dl6],  a  heterogeneous  multi¬ 
scale  method  based  on  the  discontinuous  Galerkin  method  is  developed  by  Ghen, 
E  and  Shu.  This  paper  demonstrates  the  good  potential  of  discontinuous  Galerlin 
method  in  multiscale  calculations. 
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WENO  (weighted  essentially  non-oscillatory)  finite  difference  and  finite  volume 
methods  have  been  extensively  developed  and  applied  during  this  period.  In  [a3], 
these  methods  are  compared  with  the  discontinuous  Galerkin  method  and  a  guideline 
is  given  as  to  when  each  method  has  its  unique  advantage.  In  [a5],  a  technique  is 
introduced  to  treat  the  appearance  of  negative  linear  weights  in  WENO  reconstruc¬ 
tions  while  maintaining  the  stability  of  the  approach.  In  [a8,al5,a24]  and  [b2],  WENO 
methods  are  developed  to  solve  models  in  semiconductor  device  simulations,  including 
the  direct  numerical  simulation  via  the  Boltzmann-Poisson  equations,  which  is  very 
difficult  because  of  the  high  dimensions  (a  two  dimensional  simulation  would  involve 
2  space  dimensions,  3  phase  dimensions  plus  time,  i.e.  5-1-1  dimensions).  In  [alO], 
high  order  central  WENO  schemes  are  developed  and  analyzed.  A  comprehensive 
survey  of  WENO  schemes  and  DO  methods  is  given  in  [all].  In  [al4],  a  high  order 
WENO  method  is  developed  for  solving  the  Hamilton-Jacobi  equations  on  arbitrary 
triangulations.  In  [al7],  a  multi-domain  finite  difference  WENO  method,  which  can 
be  used  on  quite  general  geometry  and  yet  is  much  less  expensive  than  finite  volume 
WENO  method,  is  developed  and  applied  to  computational  fluid  dynamics  problems. 
In  [al9,a22]  and  [b5],  the  resolution  properties  of  high  order  WENO  schemes  when 
used  on  problems  with  both  shocks  and  complicates  smooth  structures  are  explored, 
and  it  is  concluded  that  higher  order  WENO  schemes  still  have  advantages  in  obtain¬ 
ing  a  comparable  resolution  with  smaller  CPU  cost.  In  [a20],  a  Lax-Wendroff  type 
time  discretization  procedure  is  developed  for  finite  difference  WENO  schemes,  which 
could  be  more  efficient  than  the  traditional  Runge-Kutta  time  discretization  under 
certain  circumstances.  In  [a21],  a  WENO  solver  is  developed  for  a  multi-class  MWR 
traffic  model.  In  [a26],  a  WENO  scheme  is  developed  for  a  particle-fluid  two  phase 
problem.  In  [a32],  a  WENO  scheme  is  developed  for  cosmological  hydrodynamic  prob¬ 
lem  in  astrophysics.  In  [dll,dl5],  well  balanced,  high  order  finite  difference  WENO 
schemes  are  developed  which  can  maintain  exactly  certain  steady  state  solutions  and 
at  the  same  time  are  genuinely  high  order  accurate  for  general  solutions  of  the  shallow 
water  equations  and  the  hyperbolic  model  of  chemosensitive  movements. 

In  spectral  methods.  Carpenter,  Gottlieb  and  Shu  [al3]  proved  conservation  and 
convergence  of  a  class  of  global  schemes  including  the  spectral  methods  and  compact 
methods. 
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ods  and  simulation  for  particle-fluid  two  phase  flow  problems,  Computers  and 
Mathematics  with  Applications,  v47  (2004),  pp.  1437-1462. 

27.  Y.  Guo,  C.-W.  Shu  and  T.  Zhou,  The  dynamics  of  a  plane  diode,  SIAM  Journal 
on  Mathematical  Analysis,  v35  (2004),  pp. 1617-1635. 

28.  B.  Cockburn,  F.  Li  and  C.-W.  Shu,  Locally  divergence-free  discontinuous  Galerkin 
methods  for  the  Maxwell  equations.  Journal  of  Computational  Physics,  vl94 
(2004),  pp.588-610. 

29.  Q.  Zhang  and  C.-W.  Shu,  Error  estimates  to  smooth  solutions  of  Runge-Kutta 
discontinuous  Galerkin  methods  for  scalar  conservation  laws,  SIAM  Journal  on 
Numerical  Analysis,  v42  (2004),  pp. 641-666. 

30.  D.  Levy,  C.-W.  Shu  and  J.  Yan,  Local  discontinuous  Galerkin  methods  for  non¬ 
linear  dispersive  equations.  Journal  of  Computational  Physics,  vl96  (2004), 
pp.751-772. 

31.  Y.  Xu  and  C.-W.  Shu,  Local  discontinuous  Galerkin  methods  for  three  classes  of 
nonlinear  wave  equations.  Journal  of  Computational  Mathematics,  v22  (2004), 
pp.250-274. 

32.  L.-L.  Feng,  C.-W.  Shu  and  M.  Zhang,  A  hybrid  cosmological  hydrodynamic/N- 
body  code  based  on  a  weighted  essentially  non- os  dilatory  scheme,  Astrophysical 
Journal,  v612  (2004),  pp.1-13. 


(b)  Papers  published  in  non-peer-reviewed  journals  or  in  conference  proceedings 


1.  C.-W.  Shu,  A  survey  of  strong  stability  preserving  high  order  time  discretiza¬ 
tions,  in  Gollected  Lectures  on  the  Preservation  of  Stability  under  Discretization, 
D.  Estep  and  S.  Tavener,  editors,  SIAM,  2002,  pp. 51-65. 

2.  J.A.  Carrillo,  1.  Gamba,  O.  Muscato  and  C.-W.  Shu,  Gomparison  of  Monte 
Garlo  and  deterministic  simulations  of  a  silicon  diode,  in  Transport  in  Transi¬ 
tion  Regimes,  N.B.  Ben  Abdallah,  A.  Arnold,  P.  Degond,  1.  Gamba,  R.  Glassey, 
C.  Levermore  and  C.  Ringhofer,  editors,  IMA  Volumes  in  Mathematics  and  Its 
Applications,  vl35,  Springer-Verlag,  New  York,  2003,  pp. 75-84. 

3.  C.-W.  Shu,  An  overview  on  high  order  numerical  methods  for  convection  dom¬ 
inated  PDEs,  in  Hyperbolic  Problems:  Theory,  Numerics,  Applications,  T.Y. 
Hon  and  E.  Tadmor,  editors,  Springer-Verlag,  Berlin,  2003,  pp. 79-88. 
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4.  Y.-T.  Zhang  and  C.-W.  Shu,  Third  and  fourth  order  weighted  ENO  sehemes  for 
Hamilton- Jaeobi  equations  on  2D  unstruetured  meshes,  in  Hyperbolie  Problems: 
Theory,  Numeries,  Applieations,  T.Y.  Hon  and  E.  Tadmor,  editors,  Springer- 
Verlag,  Berlin,  2003,  pp. 941-950. 

5.  Y.-T.  Zhang,  J.  Shi,  C.-W.  Shu  and  Y.  Zhou,  Resolution  of  high  order  WENO 
sehemes  and  Navier-Stokes  simulation  of  the  Rayleigh- Taylor  instability  prob¬ 
lem,  in  Computational  Fluid  and  Solid  Mechanics  2003,  K.J.  Bathe,  Editor,  the 
Proceedings  of  the  Second  MIT  Conference  on  Computational  Fluid  and  Solid 
Mechanics,  June  17-20,  2003,  volume  1,  pp. 1216-1218,  Elsevier  Science. 

6.  B.  Cockburn,  F.  Li  and  C.-W.  Shu,  Diseontinuous  Galerkin  methods  for  equa¬ 
tions  with  divergenee-free  solutions:  preliminary  results,  in  Computational  Fluid 
and  Solid  Mechanics  2003,  K.J.  Bathe,  Editor,  the  Proceedings  of  the  Second 
MIT  Conference  on  Computational  Fluid  and  Solid  Mechanics,  June  17-20, 
2003,  volume  2,  pp. 1900-1902,  Elsevier  Science. 


(d)  Manuseripts  submitted,  but  not  published 


1.  F.  Li  and  C.-W.  Shu,  Loeally  divergenee-free  diseontinuous  Galerkin  methods 
for  MHD  equations.  Journal  of  Scientific  Computing,  to  appear. 

2.  M.  Zhang  and  C.-W.  Shu,  An  analysis  of  and  a  eomparison  between  the  diseon¬ 
tinuous  Galerkin  and  the  speetral  finite  volume  methods.  Computers  and  Fluids, 
to  appear. 

3.  Y.  Ha,  C.  Gardner,  A.  Gelb  and  G.-W.  Shu,  Numerieal  simulation  of  high  Maeh 
number  astrophysieal  jets  with  radiative  eooling.  Journal  of  Scientific  Gomput- 
ing,  to  appear. 

4.  W.-S.  Don,  D.  Gottlieb,  G.-W.  Shu,  O.  Schilling  and  L.  Jameson,  Numerieal 
eonvergenee  study  of  nearly-ineompressible,  inviseid  Taylor-Green  vortex  flow. 
Journal  of  Scientific  Gomputing,  to  appear. 

5.  J.  Ryan,  G.-W.  Shu  and  H.  Atkins,  Extension  of  a  post-proeessing  teehnique  for 
the  diseontinuous  Galerkin  method  for  hyperbolie  equations  with  applieation  to 
an  aeroaeoustie  problem,  SIAM  Journal  on  Scientific  Gomputing,  to  appear. 

6.  J.  Qiu  and  G.-W.  Shu,  Runge-Kutta  diseontinuous  Galerkin  method  using  WENO 
limiters,  SIAM  Journal  on  Scientific  Gomputing,  to  appear. 

7.  J.  Qiu  and  G.-W.  Shu,  Hermite  WENO  sehemes  and  their  applieation  as  lim¬ 
iters  for  Runge-Kutta  diseontinuous  Galerkin  method  II:  two  dimensional  ease, 
Gomputers  and  Fluids,  to  appear. 
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8.  J.  Qiu  and  C.-W.  Shu,  Hermite  WENO  schemes  for  Hamilton- Jacobi  equations, 
Journal  of  Computational  Physics,  to  appear. 

9.  F.  Li  and  C.-W.  Shu,  Reinterpretation  and  simplified  implementation  of  a  dis¬ 
continuous  Galerkin  method  for  Hamilton- Jacobi  equations,  Applied  Mathemat¬ 
ics  Letters,  to  appear. 

10.  J.  Qiu,  M.  Dumbser  and  C.-W.  Shu,  The  discontinuous  Galerkin  method  with 
Lax-Wendroff  type  time  discretizations,  submitted  to  Computer  Methods  in 
Applied  Mechanics  and  Engineering. 

11.  F.  Filbet  and  C.-W.  Shu,  Approximation  of  hyperbolic  models  for  chemosensitive 
movement,  submitted  to  SIAM  Journal  on  Scientific  Computing. 

12.  Y.  Xu  and  C.-W.  Shu,  Local  discontinuous  Galerkin  methods  for  the  Kuramoto- 
Sivashinsky  equations  and  the  Ito-type  coupled  KdV  equations,  submitted  to 
Computer  Methods  in  Applied  Mechanics  and  Engineering. 

13.  D.  Levy,  S.  Nayak,  C.-W.  Shu  and  Y.-T.  Zhang,  Gentral  WENO  schemes  for 
Hamilton- Jacobi  equations  on  triangular  meshes,  submitted  to  SIAM  Journal 
on  Scientific  Computing. 

14.  Y.  Xu  and  C.-W.  Shu,  Local  discontinuous  Galerkin  methods  for  nonlinear 
Schrodinger  equations,  submitted  to  Journal  of  Computational  Physics. 

15.  Y.  Xing  and  C.-W.  Shu,  High  order  finite  difference  WENO  schemes  with  the 
exact  conservation  property  for  the  shallow  water  equations,  submitted  to  Jour¬ 
nal  of  Computational  Physics. 

16.  S.  Chen,  W.  E  and  C.-W.  Shu,  The  heterogeneous  multi-scale  method  based 
on  the  discontinuous  Galerkin  method  for  hyperbolic  and  parabolic  problems, 
submitted  to  Multiscale  Modeling  and  Simulation:  A  SIAM  Interdisciplinary 
Journal. 
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